Controlled transfer of quantum amplitude via modulation of a potential barrier: 

numerical study in a model of SQUID 
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We numerically integrate the time-dependent Schrodinger equation in a single-degree-of-freedom 
model of SQUID with a variable potential barrier between the basis flux states. We flnd that linear 
superpositions of the basis states, with relatively little residual excitation, can be formed by pulsed 
modulations of the barrier, provided the pulse duration exceeds the period of small oscillations of 
the flux. Two pulses applied in sequence exhibit strong interference effects, which we propose to use 
for an experimental determination of the decoherence time in SQUIDs. 
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Ability to manipulate linear superpositions of quantum 
states is expected to lead to significant advances in infor- 
mation processing, the paradigm that became known as 
quantum computing In addition, if the basis states 
are macroscopically distinct, one may be able to gain use- 
ful insights into how quantum coherence is destroyed (or 
maintained) in the macroscopic world . 

Coherent superpositions of macroscopically distinct 
states have been obtained in two recent experiments [plH] 
on superconducting quantum interferometers (SQUIDs). 
The SQUIDs were in the regime where the basis states, 
corresponding to different values of the magnetic flux, 
were separated by a potential barrier. In the first ex- 
periment a coherent superposition was obtained by 
exciting a SQUID with a pulse of microwave radiation, 
thus bringing the system closer to the top of the poten- 
tial barrier separating the basis states. In the second 
experiment, the barrier was low by design 

For quantum information processing, on the other 
hand, it is essential that the system can be excited and 
subsequently deexcited so as to produce a controlled su- 
perposition of stable, low-lying states. The energy of 
these states, after deexcitation, should be much lower 
than the height of the potential barrier separating them, 
to ensure that no unwanted transitions occur afterwards. 
The protocols used in refs. j^,^ do not allow for such 
deexcitation. 

It is clear that a system with a fixed low-height barrier 
will not be able to satisfy the above requirement. Thus, 
it is natural to look at systems with variable potential 
barriers: lowering the barrier, by some external means, 
will correspond to excitation, while restoring it to the 



original height, to deexcitation. The combined effect of 
the two operations will be writing a stable quantum su- 
perposition to the device. 

While at present it is far from obvious that SQUIDs 
will become the underlying technology for quantum com- 
puters, these systems do allow modulation of the barrier 
height (which in this case is the Josephson coupling en- 
ergy) by an external agent, such as electric current or 
magnetic field. Indeed, the height of the barrier was ad- 
justable in the experiment of ref. 

The aim of the present work was to find, through nu- 
merical integrations, if modulation of the potential bar- 
rier in such a system, by a pulse of current or field, can 
be useful for obtaining controlled superpositions of basis 
states, and what the required durations of pulses may 
be.0 As a model, we used the reduced, single-degree- 
of-freedom, description of a SQUID with the parameters 
adopted from ref. |^ . 

Our main result is that rather clean transfers of the 
quantum amplitude in such a system are possible within 
a few tens of picoseconds. (In Sect. 3 we present a mea- 
sure of how "clean" the transfer is.) This time scale de- 
termines the switching time, to, of the qubit. The number 
of useful operations that the qubit can perform is limited 
by the ratio td/ta, where is the decoherence time, as 
well as by another factor, which we discuss in Sect. 3. 

Decoherence is an intrinsically many-body effect 
and cannot be studied in the single-degree-of-freedom 
model. That model, however, helps to identify interfer- 
ence phenomena that can be used for an experimental 
measurement of the decoherence time in SQUIDs. In 
Sect. 4, we propose to use, for that purpose, interference 
between two consecutive pulses. This method measures 



^ Our method is different from that recently proposed in 
* On leave of absence from Physics Department, Ain Shams =^ef. | , where one changes the biasing flux through a SQUID, 
University Cairo Egypt while keeping the Josephson energy fixed (and small). 
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directly the decoherence time, as opposed to the rate of 
dissipation, which can be measured by other methods |^] 
but then has to be related to td- 



II. THE MODEL 

We start with the usual model of a SQUID as a system 
with a single degree of freedom, corresponding to the 
magnetic flux through the loop. Fluxes are conveniently 
measured in units of the flux quantum <l>o. If we denote 
the total flux in these units as y, and the biasing flux by 
y, the Schrodinger equation for the wave function '^{y, t) 
will read 



-^dl + EL{y - yf - Ej{t) cos(2^y) 



(1) 



where Eq is the charging energy due to the junction's 
capacitance C: Ec — e^/2C; El is the magnetic field 
energy, due to the inductance L: El = ^^/2L; and Ej 
is the Josephson coupling energy. In our case, the latter 
depends on time, causing a variation in the height of the 
potential barrier separating states with different values 
of the flux. 

In this paper, we only present results for the case when 
the SQUID is biased by exactly half of the flux quantum: 



y=l/2 



(2) 



although results for an arbitrary bias can be obtained 
similarly. In the case (|^), the potential in eq. is 
symmetric about the point y — 1/2; in particular, it 
has two degenerate minima. Accordingly, we introduce a 
shifted variable: 



x^y-1/2. 
We also define a new time variable t via 
t = hT . 



(3) 



(4) 



Frequencies, with respect to this new time, are measured 
in energy units, for which we use degrees Kelvin. There- 
fore one unit of time r. At = 1, corresponds to 



At: 



n 



fcs X 1 K 



7.64 ps 



(5) 



of the physical time. Finally, we shift the potential down 
by a constant equal to the unperturbed Josephson energy 
Eq. This corresponds to using the wave function tp = 
\I'exp(ii?ot/?i). Eq. ^) then leads to 



-^dl + Elx^ + Ej{t) cos(27rx) 



Ef) 




FIG. 1. Grayscale plot of the probability Pl for the system 
to remain in the left well after a Gaussian pulse of magnitude 
A and duration tq. Darker color means larger Pl, i.e. smaller 
transfer of the amplitude to the right. 

We use parameters of the SQUID described in ref. 
@: Ec = Q mK, El = 645 K, Eq = 76 K. The time- 
dependent Josephson energy is taken in the form 



Ej{T)/Ea = 1 - Aexp[-(r - /t^] 



(7) 



This corresponds to a Gaussian (in time) pulse of current 
through the circuit that controls the Josephson coupling. 
If the time Tc, corresponding to the center of the pulse, is 
noticeably larger than the duration of the pulse, tq, the 
initial value of Ej, Ej{0), is practically indistinguishable 
from Eq. This initial value corresponds to a potential 
barrier in high enough for the tunneling between dif- 
ferent flux states to be inefficient. Accordingly, we can 
talk about the left and right ground states, which are 
the lowest states in each well. Near t = Tc, however, the 
barrier is lower, and if 



A>A„^l- El/2tt^Eo = 0.57 



(8) 



(6) 



it disappears altogether. At large times, the barrier 
comes back to its original height. 



III. NUMERICAL RESULTS 

We numerically integrated eq. (^, with a Gaussian 
pulse of the form (^. In the preparation of each run, we 
started by relaxing the system to its true ground state 
wavefunction (which is a superposition of the left and 
right ground states). Then we selected (and normalized) 
the left half of that wavefunction as our initial state. Af- 
ter that we numerically followed the system's evolution. 

The magnitude of the pulse. A, and the duration tq 
were varied, with the goal to achieve either a complete 
transfer of the amplitude between the left and right wells, 
or a superposition state. We concentrated primarily on 
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values of tq large enough for the transition to be adia- 
batic, so that in the end the system was not significantly 
excited beyond the left and right ground states. 

Throughout the course of the evolution, we monitored 
a number of quantities. One of these was the probability 
Pl to find the system in the left well after the pulse was 
completed. It is shown in Fig. |l| for a range of values of A 
and Tq. For values of A larger than the critical value (H), 
we observe that Pl is quasiperiodic with respect to tq, 
and the period decreases as A increases. We interpret this 
by noting that for such values of A, the potential during 
the pulse has a minimum at a; = 0, about which the 
system can oscillate. Duration of the pulse will determine 
in what phase of these oscillations the system will be 
deexcited. 

Theoretically, we expect oscillations of Pl , due to tun- 
neling back and forth between the two wells, even when A 
is well below the critical value. However, the correspond- 
ing timescale is too long for such oscillations to show up 
in the figure. 

Fig. |l| indicates that multiple choices of the parame- 
ters may lead to the same final Pl ■ The presence of deco- 
herence in any realistic system suggests that one should 
choose the smallest possible tq, so as to achieve the max- 
imal number of quantum memory switches before the de- 
vice decoheres. On the other hand, a too small tq will 
result in overexcitation, and after many such switches 
the control of the system will be lost. (One may con- 
template reducing excitation through cooling; however, 
any cooling apparatus is likely to become a significant 
source of decoherence.) Thus, in addition to the quality 
factor Q = td/to, where td is the decoherence time, we 
introduce a "fidelity factor" 



F^iEb^ Eg)/AE , 



(9) 



where Ef, is the energy at the top of the unperturbed 
potential barrier. Eg is the ground state energy, and AE 
is the increase in energy after one switch. Decreasing to to 
achieve a larger Q leads, via the time-energy uncertainty, 
to an increase in AE and a smaller F, not unlike how 
increasing the clock rate in ordinary computers leads to 
overheating. 

To study the "fidelity" of transitions in more detail, we 
now focus on two values of A: A = 0.53 and 0.59 (corre- 
sponding to two sections of Fig. 1). Both are close to the 
threshold (H). However, the first is below that threshold, 
so tunneling effects are present. In Fig. ^ we show the 
corresponding probabilities Pl together with the rescaled 
final energy E of the system. In our case, Et, — (be- 
cause of the way the potential has been shifted), and Eg 
is negative {Eg — —41.1), so the plotted quantity E/\Eg\ 
is simply related to the fidelity parameter M): 



E 1 



E 



(10) 



For example, for A = 0.59 and tq = 5, we obtain Pl = 
0.004 (so the amplitude is almost completely transferred 



to the right) and E = —39.8, corresponding to F 30. 
Increasing tq to tq = 35 {Pl = 0.09) results in F ^ 100. 
We also observe that, perhaps contrary to intuition, for 
short pulses the fidelity is higher when the magnitude 
is slightly above the threshold (|8|), rather than slightly 
below. 
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FIG. 2. Probability Pl, as a function of the pulse duration, 
for A — 0.59 (solid line) and A = 0.53 (long dashes). The ra- 
tio of the final state energy to the magnitude of the ground 
state (initial) energy is shown by short dashes and dots, re- 
spectively. This ratio is a measure of the excitation energy 
supplied to the system by the pulse. 

Finally, Fig. ^ illustrates how the transition develops 
in time. In this case, it is predominantly an overbarrier 
transition, rather than tunneling. In the full many-body 
theory of the SQUID, such overbarrier motion of the flux 
will be associated with dissipation and decoherence. It 
remains to see if the time required for the transition can 
be made sufficiently short for these effects not to present 
a problem. 
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FIG. 3. A few profiles of the probability density, showing 
how the transition occurs for A = 0.59 and tq = 5.1. The 
maximum of \ip\'^ moves from left to right as time increases. 
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IV. TWO-PULSE INTERFEROMETRY 



pm = \rLoum\' + \rLium\' 



We now discuss the behavior of the system under two 
consecutive pulses and the possibihty to use such two- 
pulse sequences for an experimental determination of the 
decoherence time. Here we discuss decoherence after a 
transition, as opposed to decoherence during it, men- 
tioned above. In its general outline, the method is simi- 
lar to the two-pulse method used in studies of magnetic 
resonance [Q: in either case, the final state depends on 
the phase accumulated by some part of the wavefunction 
during the interval between the pulses. We now describe 
our method in more detail. 

Suppose that the system is originally in its left ground 
state, which we call ipLo- The first pulse will transfer part 
of the amplitude to the right well. Because the transi- 
tion inevitably introduces some excitation, the state on 
the right will not be the precise right ground state, ipROj 
and the state on the left will not remain the precise left 
ground state. For simplicity, and during this discussion 
only, let us assume that only the admixture of the first 
excited states, ipLi and tpm, is substantial, while higher 
excited states can be neglected. Then, the state after the 
first pulse is 

1p = CLO^LO + CLl-ipLl + CrqiI^ro + CRilpRi . (11) 

Time evolution of this state is given by 

Ipit) = [CLO^LO + CRolpBO + 

(cLiV'Li + CRiil^m) exp(-iLji)] eiip{~iEgt/h) , (12) 

where u is the angular frequency for transitions between 
the ground and the first excited states. (The potential is 
assumed left-right symmetric.) 




FIG. 4. Probability that the system remains in the left well 
after two consecutive pulses separated by At. The pulses were 
identical with magnitude A — 0.59 and duration tq = 11.35. 

If at some time t we apply a second pulse, whose action 
on the state (12) is represented by some unitary operator 
U, and then measure the probability to find the system 
somewhere in the left well, the result will be 



which in general is an oscillating function of t with fre- 
quency LO. 

The amplitude of such oscillations can be quite large, 
as we show in Fig. IJ. This figure is based on numerical 
integration and does not involve the simplifying assump- 
tion (0). 

In practice, the system is coupled to the environment, 
and this coupling introduces a certain amount of decoher- 
ence. Suppose that evolution of the environment depends 
only on which well the system is in. We then consider 
two histories of the environment, xl and XRj and instead 
of eq. ( p^ ) obtain 

= {cLoi^LOXL(t) + CRotpRQXR{t) 

+ [cLiipLiXL{t) + CRi^RiXR{t)] exp(-iwt)} 
X exp{—iEgt/h) . 

The probability to observe the system in the left well 
after the second pulse is now 

p'Lit) = E [irLox^umf + \rLix*numf] , 

n 

where Xn is a complete system of states of the environ- 
ment. The oscillating dependence of on time is now 
modulated by the coherence factor |X/j(Oxl(OI^: which 
decreases in time. Note that even when this coherence 
factor becomes essentially zero, the oscillations of re- 
main, due to nonzero matrix elements, such as tP^qUiPli, 
between the states in the same well. 

Thus, due to decoherence, the oscillations of P^, such 
as those shown in Fig. ^, will acquire an envelope. By 
fitting the envelope with an exponential 



f{At) = ai + a2 exp(-At/td) 



(13) 



(or, perhaps, some other function that may emerge on 
theoretical grounds), one may be able to extract the de- 
coherence time td- 
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